A Mesoscopic Approach to the "Negative" Viscosity Effect in Ferrofluids 
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We present a mesoscopic approach to analyze the dynamics of a single magnetic dipole under 
the influence of an oscillating magnetic field, based on the formulation of a Fokker-Planck equation. 
The dissipated power and the viscosity of a suspension of such magnetic dipoles are calculated from 
non-equilibrium thermodynamics of magnetized systems. By means of this method we have found a 
non- monotonous behaviour of the viscosity as a function of the frequency of the field which has been 
^\ , referred to as the "negative" viscosity effect. Moreover, we have shown that the viscosity depends 

0^ ' on the vorticity field thus exhibiting non-Newtonian behaviour. Our analysis is complemented with 

0^ I numerical simulations which reproduce the behaviour of the viscosity we have found and extend the 

■ scope of our analytical approach to higher values of the magnetic field. 

c3 ■ Pacs: 75.50.Mn, 66.20.-fd, 05.40. +j 

Si 

m ■ 

! I. INTRODUCTION 
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l/^ • It is well-knovifn that the dynamics of a suspension of dipolar particles is strongly influenced by the presence of 
' an external field. Concerning the viscosity of the suspension, its behaviour as a function of the magnetic field is 
monotonous in the case of a constant field, reaching a saturation limit where the magnetic moments are oriented 
CO ' along the field. In contrast, when the field oscillates that behaviour is considerably modified to the extent that the 
contribution of the oscillating field to the effective viscosity of the suspension may become negative for frequencies of 
the field larger than the local vorticity. This phenomenon observed experimentally in and reported in refs. ||^ and 
has been referred to as the "negative" viscosity effect in ferrofluids. It reveals the presence of two regimes, one 
"j^ ■ essentially dissipative in which the variation of the viscosity is positive and other in which the energy of the oscillating 
[ field is practically transformed into kinetic energy of the particle. It is in this last regime where the variations of the 
I . viscosity become negative. 
''O ' Our purpose in this paper is to present an explanation of that effect based on a Fokker-Planck dynamics describing 
Ch , the time evolution of the probability density for the orientation of the dipolar particle. The starting point is the 
O ■ formulation of that equation and its perturbative solution. In this way we compute the different components of the 
susceptibility. The dissipated power and the viscosity follow from an analysis based on non-equilibrium thermody- 
namics. 

The general methodology we introduce accounts for the results of ref. valid in the limit of small field, and agree 
with the experimental data of ref. . Our theoretical results are also compared with numerical simulations we have 
performed, indicating that the qualitative behavior of the system is satisfactorily reproduced even for a moderate 
intensity of the oscillating magnetic field. 

The paper is organized as follows. Using linear response theory, we compute in Section II the generalized sus- 
ceptibility associated with the orientation vector of the magnetic dipoles in the ferrofluid, when this is under the 
influence of an oscillating magnetic field. Section III is devoted to obtain the contributions of the oscillating field to 
the dissipation of energy and to the effective viscosity. In Section IV we report results on the viscosity obtained from 
numerical simulations, whereas in Section V we present our main conclusions. 
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II. RESPONSE OF A MAGNETIC DIPOLE TO AN OSCILLATING FIELD 

We consider a dilute colloidal suspension of ferromagnetic dipolar spherical particles with magnetic moment, 

m = nisR, where R is an unit vector accounting for the orientation of the dipole. Each dipole is under the influence 
of a vortex flow with vorticity 51 = '^^oz, with z being the unit vector along the z-axis, and of an oscillating field 
H = He^^'^^x, with x being the unit vector along the x-axis. For t ^ Tr, the motion of the particle is overdamped 
with Tr = I /^r being the inertial time scale. Here / is the moment of inertia, = Sirrjoa^ is the rotational friction 
coefficient, with rjo the solvent viscosity, and a the radius of the particle. This time scale defines a cut-off frequency 
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LOr = Tj, ^ , such that the condition for overdamped motion is equivalent to w <IC ujr- In this case, the balance of the 
magnetic and hydrodynamic torques acting on each particle 

m X + ( ill - flp ) = 0, (1) 



together with the rigid rotor evolution equation 



^ = npxk (2) 



lead to the dynamic equation for R 

^ = {cjo^ + A(i) (^xf)} xl (3) 

Here A(t) = {msH/^r)^^^'^^ , with dp being the angular velocity of the particle. 

With allowance for Brownian motion, the stochastic dynamics corresponding to eq. is given by the Fokker-Planck 

equation, involving the probability density ^(i?, 

= (/:o + A(t)£i)4'(i?,i), (4) 

where Cq and Ci are operators defined by 

Co = — cjo 2 ■ 7? + DrTZ^ , 

(5) 

£i = 2^-f - {Rxlc)-n, 

with Dr ~ ksT/^r being the rotational diffusion coefficient, and TZ — Rx d/dR the rotational operator. Notice that 
the first and second terms on the right hand side of eq. (|^)i, correspond to convective and diffusive term, respectively. 
Moreover eq. which, according to eq. (jsj), rules the Brownian dynamics in the case of overdamped motion, is valid 
in the diffusion regime. This regime is also characterized by the condition t 3> t^, or equivalently lo <C tOr, which 
implicitly involves the white noise assumption. 

To solve the Fokker-Planck equation we will assume that Ao = |A(t)| constitutes a small parameter such that 
this equation can be solved perturbatively. Thus up to first order in A, the solution of the Fokker-Planck equation 
(i is 

= e(*-*°)^°*"(to) + / dt' X{t')e^'-''^^°Ci'^^{t') (6) 

J to 

Here ~ e^* ~*o)£ov|/0('^ _ i^'^ jg i^j^g ^ero order solution at time t' , and 

^\R,t = to) = 5{R-Ro), (7) 

with i?o being an arbitrary initial orientation. As follows from eq. the unperturbed operator Cq is composed of 

the operators TZz and 7?.^, which are proportional to the orbital angular momentum operators of quantum mechanics 
Lz and L^, respectively, and, therefore, their eigenfunctions are the spherical harmonics |^ 

nzYlm{R) = imYlr,^{R), 

n^Yi^{R)^-l{l + l)Yu,,{R). (8) 

Given that we know how TZ acts on the spherical harmonics, it is convenient to expand the initial condition in 
series of these functions, since the spherical harmonics constitute a complete set of functions which are a basis in the 
Hilbert space of the integrable functions over the unit sphere H 
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^\li,to) = S{R-ko) = E E ^™(^o)ilm(i?) • (9) 

1=0 m=-l 

Using this expansion in eq. (^), for the first order correction to the probabihty density, A^* = — ^P", we obtain 

<Xj I pt ^ 

Av[/(i?,t) di'A(i')i"/:„(^o)e(*-*')^°£ie(*'-*«)^°rz,„(i?). (10) 

Notice that the integral of A'i'{R, t) over the entire solid angle is zero, in agreement with the fact that the unperturbed 

solution ^'"(i?, is normalized. 

Since, we are interested in the asymptotic behavior we will set tp —>■ —oo. In this limit, eq. (^) becomes 

(11) 



where now 



and 



dt' X{t')e'-*-''^^° 



(12) 



(13) 



is the uniform distribution function in the unit sphere. 

^From eq. (na) the contribution of the AC field to the mean value of the orientation vector R can be obtained as 



R{t) 



1 

47r 



dt'X{t') / dRRe^*-''^^"2R-^. 



This equation can be written in the more compact form 



where the response function ||8| has been defined as 

1 



i?,(t) = / dt'Xit')xr{t~t'), 



dRR.e''^" 2R-i, 



for T > 0. 

By causality, we can write f ^ oo in the upper limit of the integral in eq. (|l5|); hence, this equation becomes 

R,{t) - X;{^)X{t). 

where Xi(w) is the generalized susceptibility, which is the Fourier transform of Xii^) 

dre"^^ I dRR,e^^°2R- oc 



From this equation we obtain the components of the susceptibility 



(14) 
(15) 
(16) 

(17) 
(18) 



2Dr 



2Dr. 



(19) 



(20) 
(21) 

The quantities Xx, and Xy^ have poles aX uj — icjo ± 2Dri. The inverse of the imaginary part of these poles, {2Dr)~^, 
defines the Brownian relaxation time. 
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III. NON-EQUILIBRIUM THERMODYNAMICS OF THE RELAXATION PROCESS 



Our purpose in this section is to compute the energy dissipated during the relaxation process of the magnetization 
and the rotational viscosity of the suspension. The starting point is the entropy production corresponding to the 
relaxation of the magnetization, as given in ref. |^ 

Here M is the magnetization of the suspension and H the magnetic field. Moreover, Heq is the magnetic field related 
to the instantaneous value of AI. The primes indicate that the corresponding quantities have been computed in the 
frame of reference rotating with the fluid. 

The entropy production can alternatively be written in terms of the corresponding quantities in the laboratory 
frame. Using the relation between the temporal derivatives of the magnetization in both frames 

dM' dM 1-; 

— — = — -17 X M . (23) 

dt dt 2 ^ ' 



from eq. ([2^) one obtains 

• (Heg -H) + ^nxM- [H,, - H) 

where now the different quantities refer to the laboratory system. Notice that in our case M — c mgR, where c is 
the concentration of particles. Moreover, eq. (|2^) written in this way ensures the frame material invariance of the 
entropy production. 

The linear law inferred from this expression coincides with the relaxation equation postulated by Shliomis 
provided we identify the phenomenological coefficient with the inverse of the Brownian relaxation time. 
The form of eq. ( p^ ) suggests the following decomposition 

a = (7d + (Jv ■ (25) 

The first contribution 



1 dM ,^^^ldM^ 



accounts for the entropy production coming from Debije relaxation, whereas the second 

av = X M) ■ (Heg -H) = • {M x H) (27) 

is related to the viscous dissipation. To obtain these expressions use has been made of the relation M = KHeq, with 
K being the static susceptibility, and of the fact that M remains constant which implies [dM /dt) ■ Heq — 0, and 

(nxM)- Heq = 0. 

The power dissipated in a period of the field follows from the entropy production we have computed. For the 
different contributions one has 



271- Jo 



27r/w 

Taadt , {a = D,V) . (28) 



For the Debije contribution we obtain 



PdH - 6?7o0^ IrnXxH Mt? = (29) 
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where we have used eqs . ( p7| ) and (|T9|), whith cf) = (47r/3)a^c being the volume fraction of particles. Similarly, by 
using eqs. and ( |20| ) the viscous part yields 



(30) 



rio(j>u!oXo 



LOq — LU 



This last contribution introduces the rotational viscosity defined through the relation 
Combining eqs. (^) and (^Tj) one infers the value of the rotational viscosity 



Vr 



-T]o(j)UJo ^Ag 



UJq — LU 



(31) 



(32) 



which represents the contribution of the rotational degrees of freedom of the dipoles to the effective viscosity of the 
suspension given by 



Veff =V + Vr 



(33) 



where r] is given by Einstein law = 770(1 + 5 /2(f)), [0. In view of the former results, the total dissipation P{lu) = 
Pd(w) +Pv(w) is 



(wo + ujY 



(34) 



Notice that although the contributions to the dissipation Pd{'-^) and Py(aj) may achieve negative values, the total 
dissipation remains always positive, in accordance with the second law. 

The linear law derived from the entropy production (24) is valid for situations closed to equilibrium, when the 
distribution function is not too different from the equilibrium distribution. For larger deviations, this approach and 
equivalently the relaxation equation of Shliomis are no longer valid. One then should employ the Fokker Planck 
description we have introduced in section II. In this way, results for higher values of the field could also be obtained 
by means of the for malism we have developed. 

IV. NUMERICAL SIMULATIONS 

In order to check the validity of our results and explore the behaviour of the system for higher values of the 
oscillating field we have performed numerical simulations by using a standard second-order Runge-Kutta method for 
stochastic differential equations [T^JT^ . To this purpose we have considered the Langevin equation corresponding to 
eq. (i: 



dR 

dt 



X i?. 



(35) 



where ^(t) is Gaussian white noise with zero mean and correlation {^i{t)^j{t + r)) = 2DrSijS{T). 

^From the previous equation one can easily compute the mean angular velocity of the particles by averaging over 
several realizations of the noise. This quantity then gives the rotational viscosity 



|,.*(i-2' 

2 V Wo 



(36) 



To obtain this expression we will first rewrite the entropy production in terms of P", the axial vector related to the 
antisymmetric part of the pressure tensor 



1 

T' 



p" -n 



(37) 
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from which one derives the phenomenological law 



P° ^ rjr^ . (38) 
The balance of torque densities, which is achieved in the asymptotic regime |]lo|| , 

2P" = M X i? = 6r]o(l){^n - dp) (39) 

together with eqs. (^gf), and ( ^ ) then leads to expression (^6|). 

In Figs. (|l]) and (^ we show the results obtained for the normalized rotational viscosity, rjr /riQ(j)\^, as a function 
of the frequency of the applied field and as a function of the angular velocity of the fluid, respectively, for different 
intensities of the oscillating field. These figures clearly illustrate how for low intensities of the oscillating field the 
results obtained through numerical simulations can accurately be reproduced by linear response theory, as well as for 
higher values, linear response theory is still qualitatively correct. The crossover point from positive to negative values 
of the rotational viscosity seems not to depend significantly on the amplitude of the oscillating field. It is interesting 
to realize that, at high frequencies all the simulation curves match the linear response theory curve. Thus, in this 
frequency regime linear response theory provides an accurate description of the phenomenon. Fig. (^ also makes the 
non-Newtonian character of the fiuid manifest. Notice that these effects are more pronounced for angular velocities 
of the fluid near the angular frequency of the applied field, i.e., for high and low values of the vorticity the rotational 
viscosity goes to zero and a constant negative value, respectively, whereas for intermediate values it depends on the 
particular value of the vorticity. From both figures, one concludes that the vorticity plays just the opposite role to 
the frequency of the field. 



V. CONCLUSIONS 



In this paper we have presented a mesoscopic approach to explain the "negative" viscosity effect ocurring in 
suspensions of magnetic particles. The dynamics of the magnetic moment in the oscillating field is described by 
means of a Fokker-Planck equation which can be solved perturbatively. This equation gives rise to a hierarchy of 
equations fl^ for the different moments, describing the relaxation of the magnetic moment. The first equation of the 
hierarchy, when linearized in the field , agrees with the phenomenological equation obtained from non-equilibrium 
thermodynamics and the corresponding one postulated by Shliomis. 

Following this procedure, we have been able to compute the dissipated power and the viscosity which is a non- 
monotonous function of the frequency of the field. Up to first order in the field our results agree with the corresponding 
ones of ref. based on the phenomenological equation proposed by Shliomis. 

The phenomenological approach, dealing with the phenomenological law derived from the entropy production, is 
strictely valid at small fields and is no longer correct when larger deviations from the equilibrium distribution occur. 
If we are interested in the response of the system to larger values of the field the correct approach is the one based 
upon the Fokker-Planck equation we have proposed in section II. This approach involves the evolution equation for 
higher order moments of the probability distribution in the hierarchy. After introducing a decoupling approximation, 
the solution for the first moment can be used to compute the energy dissipation that allows one to define an effective 
viscosity. 

In order to go beyond the linear regime we have performed numerical simulations. This numerical analysis enables 
us to discern about the validity of the linear response theory treatment. Our conclusion is that qualitatively linear 
response theory may provide a reasonably explanation of the phenomenon. 
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FIG. 1. Normalized rotational viscosity, rir/rio4>\^, as a function of the frequency of the applied field for 
Ao = {12.8, 6.4, 3.2, 1.6, 0.8, 0.4} and linear response theory (LRT). The values of the remaining parameters are Dr = 0.5 
and Wo = 27r. 
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FIG. 2. NormaJized rotational viscosity, rjr/riocpXo, as a function of the angular velocity of the fluid for 

Ao = {12.8, 6.4, 3.2, 1.6, 0.8, 0.4} and linear response theory (LRT). The values of the remaining parameters are Dr = 0.5 
and w = 277. 
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